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(54) Method and apparatus for measuring power output and for measuring tension and 
vibrational frequency of a elongate flexible member 



(57) The invention broadly provides a method and 
apparatus for measuring the power output of one pow- 
ering a vehicle or other system driven by the elongate 
flexible member such as chain. In a preferred embodi- 
ment described, the vehicle is a bicycle. The method of 



the invention comprises measuring the speed and ten- 
sion of the drive chain and calculating the power output 
therefrom. The apparatus of the invention includes a 
chain speed sensor, a chain tension sensor, and elec- 
tronic processing means to calculate and display the 
power output based upon the chain measurements. 
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Description 
Field of the Invention 

5 [0001] This invention relates generally to a method and apparatus for measurement of power output of one powering 
a system driven by an elongate flexible member such as chain. The system preferably is a bicycle. More specifically 
it relates to an apparatus and method for the measurement of power output of a bicyclist. 

Background of the Invention 

10 

[0002] It is often desirable to measure power output of a bicyclist or a rider of other human powered vehicle. Such 
information is useful both for the casual rider as well as the high performance trainer. A bicyclist can use this information, 
among other things, to monitor increases in cardiovascular fitness, quantify maximum output and power, optimize 
pedaling cadence, and determine maximum efficiency. Such measurements are also useful to those seeking to reduce 

15 body fat though exercise or to replenish food during extended rides. While power measuring devices are available on 
stationary exercise equipment commonly found in health clubs and the like, the size of currently available power meas- 
uring devices makes it impractical to install them on actual human powered vehicles. Some force and power meas- 
urement equipments use strain gage sensors installed on the crank or the pedal of a bicycle. This attaching typically 
by gluing must be done carefully with special instructions. The sensor also needs calibration after attaching. One 

20 presented way to measure power is also done by sensoring the shoe of a bicyclist. This is not, however, practical in 
every day use. 

[0003] In addition to human powered vehicles, power measurements would also be useful for motor driven devices. 
For example, in motor assisted bicycles, chain tension and/or human power output can be used to determine the 
relative instantaneous output of the motor. For motor-driven devices such as motorcycles, chain tension and power 
25 output can be used to monitor the load on the engine, which relates information on fuel consumption, mechanical wear, 
and efficiencies. 

[0004] What is needed then, is a method and device which is compact in size and operatively arranged to measure 
power output of a bicyclist. Other quite necessary needs are that sensor unit can be easily attached and signal to be 
measured is easy to find with a good signal to noise ratio. On the other hand a measurement device could be taken 
30 to use without a difficult calibration procedure. 

Summary of the Invention 

[0005] The invention broadly provides a method and apparatus for measuring the power output of one powering a 
35 bicycle or other system driven by an elongate flexible member. In a preferred embodiment described, the vehicle is a 
bicycle. 

[0006] The method of the invention comprises measuring the speed and vibration of the elongate flexible member 
and calculating the power output therefrom. The apparatus of the invention includes a speed sensor, a vibration sensor, 
and electronic processing means to calculate and display the power output based upon the measurements made 

40 relating to elongate flexible member. 

[0007] An apparatus for measuring the power transmitted by an elongate flexible member in this invention comprises 
vibration sensor that measures the vibration frequency of the elongate flexible member, a sensor that directly or indi- 
rectly measures the speed of the elongate flexible member; and a computer that determines the power transmitted by 
the elongate flexible member by receiving the vibration measurement and the speed measurement directly or indirectly 

45 measured. 

[0008] A method for measuring the power transmitted by an elongate flexible member in this invention comprises 
the steps of measuring the vibration frequency of the elongate flexible member, directly or indirectly measuring the 
speed of the elongate flexible member, and receiving the vibration measurement and the directly or indirectly measured 
speed measurement in a computer, and finally generating therefrom a signal representing the power transmitted by 

so the elongate flexible member. In another embodiment, a method for determining the tension in a elongate flexible 
member on a vehicle such as a bicycle driven by the elongate flexible member comprises the steps of identifying a 
vibrating section of the elongate flexible member, providing a vibration sensor on the vehicle in proximity to the vibrating 
section of the elongate flexible member, measuring the frequency of vibration of the vibrating section of the elongate 
flexible member by means of the sensor during operation of the vehicle, and utilizing the frequency of vibration and 

55 the known length of the vibrating section of the elongate flexible member to determine the tension in the elongate 
flexible member . 

[0009] In another embodiment a method for determining the vibrational frequency of a elongate flexible member on 
a vehicle such as a bicycle driven by the elongate flexible member comprises providing the vibration signal of the 
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elongate flexible member by means of a vibration sensor in proximity to the elongate flexible member, performing a 
Fast Fourier Transform on the vibration signal of the elongate flexible member, and heuristically selecting the vibration 
frequency of the elongate flexible member. 

[0010] In another embodiment an apparatus for measuring the tension in a elongate flexible member on a vehicle 

s such as a vehicle driven by the elongate flexible member comprises a vibration sensor mounted on the vehicle in 
proximity to a vibrating section of the elongate flexible member for measuring the frequency of vibration of the vibrating 
section of the elongate flexible member during operation of the vehicle, signal processing means on the vehicle oper- 
atively connected to the vibration sensor for selecting a measured frequency of vibration, and computer means on the 
vehicle operatively connected to the signal processing means for computing the tension in the elongate flexible member. 

io [0011] In another embodiment, a method for measuring the power output of a person operating a vehicle such as a 
bicycle driven by means including a elongate flexible member comprises the steps of measuring the vibration of the 
elongate flexible member by means of a sensor on the vehicle, measuring directly or indirectly the speed of the elongate 
flexible member by means of a sensor on the vehicle, computing the power transmitted by the elongate flexible member 
by means of a computer on the vehicle using the directly or indirectly measured speed of the elongate flexible member 

15 and the vibration of the elongate flexible member. 

[0012] In another embodiment, apparatus for measuring the power output of a person operating a vehicle such as 
a bicycle driven by means including a elongate flexible member comprises a vibration sensor mounted on the vehicle 
in proximity to a vibrating section of the elongate flexible member for measuring the vibration frequency of the elongate 
flexible member, a sensor mounted on the vehicle in proximity to the elongate flexible member for directly or indirectly 

20 measuring the speed of the elongate flexible member; and a computer on the vehicle and operatively connected to 
the vibration sensor and to the sensor which directly or indirectly measures the speed, said computer computing the 
power transmitted by the elongate flexible member using the measured vibration and directly or indirectly measured 
speed of the elongate flexible member. 

[0013] In another embodiment, a method for measuring the power transmitted by an elongate flexible member the 
25 method comprising the steps of measuring directly or indirectly the speed variation of the elongate flexible member, 
and receiving the direct or indirect speed variation measurement in a computer, generating a signal representing the 
power transmitted by the elongate flexible member In this embodiment vibration measurement is not been but speed 
variation measurement is used. 

[0014] In another embodiment similar to above, an apparatus for measuring the power transmitted by an elongate 
30 flexible member comprises a system having a first engagement means separated from a second engagement means, 
the elongate flexible member coupling the first engagement means with the second engagement means so that when 
the first engagement means is rotated by a power source the elongate flexible member moves in the direction of the 
rotation of the first engagement means, which in turn causes the moving elongate flexible member to rotate the second 
engagement means, sensor that directly or indirectly measures the speed variation of the elongate flexible member; 
35 and a computer that determines the power transmitted by the elongate flexible member by receiving the vibration 
measurement and the speed variation measurement directly or indirectly measured. In this embodiment vibration meas- 
urement is not been but speed variation measurement is used. 

[0015] The primary object of the invention is to provide a method and apparatus for measuring the power output of 
one powering a chain driven vehicle. The advantages of this invention are simple easy to install sensor units avoiding 
to complex calibration process and however accurate and reliable force and power measurement. Another object of the 
invention is to display power output of a bicyclist in real time as the cyclist is riding. 

[0016] These and other objects, advantages and features of the invention will become obvious to those having or- 
dinary skill in the art upon a reading of the following detailed description, read in conjunction with the drawings. 

45 Brief Description of the Drawings 

[0017] 

Figure 1 is a side view of a typical multi-speed bicycle with which the present invention is used; 
50 Figure 2 is a fragmentary side view of the rear derailleur of the bicycle shown in Figure I and a speed sensor 

mounted to the derailleur; 

Figure 3 is a cross-sectional view of the derailleur/sensor combination of Figure 2, taken generally along line 3-3 
of Figure 2; 

Figure 4 is a fragmentary top view of the vibration sensor of the invention; 
55 Figure 5 is a side view of the sensor shown in Figure 4; 

Figure 6 is a cross-sectional view taken generally along line 6-6 in Figure 4; 

Figure 7 is a schematic block diagram of the electrical and electronic circuit of the apparatus of the invention; 
Figure 8 is an electrical schematic diagram of the connectors used to connect the peripheral interface circuit to 
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the computer processing unit; 

Figure 9 is an electrical schematic diagram of part of the peripheral interface circuit of the invention; 
Figure 10 is an electrical schematic diagram of part of the peripheral interface circuit of the invention; 
Figure 11 is an electrical schematic diagram of part of the peripheral interface circuit of the invention; and, 
5 Figure 12 is a flow chart of the computer program operatively arranged to run the digital signal processor of the 

invention. 

Detailed Description of the Preferred Embodiment 

10 [0018] At the outset, it should be understood that the description which follows is directed to a preferred embodiment 
of the invention, but that the appended claims are not so limited. It should be also appreciated that like reference 
numerals on different drawing figures refer to identical structural elements. 

[0019] Figure 1 illustrates a typical multi-speed bicycle 10 which is well known in the art. The bicycle, which, in this 
embodiment is a 10-speed bicycle, comprises a frame 11 , which frame comprises top tube 34, seat tube 32, seat stay 
'5 31, chain stay 23 and down tube 33. The bicycle further comprises front fork 17, front wheel 12, rear wheel 13, seat 
14, handlebars 15, drive sprocket 16, pedal 25, chain 18, front derailleur 21, sprocket cluster 20, rear derailleur 19, 
front derailleur shift lever 28 and rear derailleur shift lever 26. 

[0020] As is well known in the art, to propel the bicycle, the rider drives the pedal, turning the front sprocket and 
causing the chain to engage and drive one of the rear sprockets. The rear sprockets, in turn, are fixedly secured to the 

20 rear axle of the rear wheel, and thereby drive the wheel to cause the bicycle to move. The front derailleur is used to 
move the chain from one of a plurality of front drive sprockets. In a 1 0-speed bicycle as shown, there are typically two 
front sprockets, each with a different number of teeth (e.g., 22 teeth on the small sprocket and 54 teeth of the large 
sprocket). Similarly, the rear derailleur is used to move the chain from one rear sprocket to another. Typically, in a 10 
speed bicycle as shown, there are five rear sprockets, each containing a different number of teeth (e.g., 12 teeth on 

25 the smallest diameter sprocket and 28 teeth on the largest diameter sprocket). 

[0021] The power output of the cyclist is equal to the amount of work done by the cyclist per unit time. One meas- 
urement of this power is given by the equation: 

30 P(watts) = chain speed (meters/sec) * chain tension (Newtons) 

[0022] Thus, it is seen that power output is calculated from measurements of chain speed and chain tension. 
[0023] It is noticed here that power output is not typically constant during a cadence cycle of rotating cranks. Power 
fluctuates due to different force and moment achieved against pedals and cranks. The minimum power or force is 
35 typically detected when one pedal is at highest position and another pedal is at lowest position. The maximum power 
or force is typically achieved when the pedals are at the same horizontal level. Because power and force is transmitted 
directly to the chain through the front sprocket, the power variation also causes chain speed and chain tension to 
fluctuate during a cadence cycle respectively. 

[0024] It will be appreciated by those skilled in the art that the present invention can be used to measure the power 
40 transmitted by any elongate flexible member, for example, a bicycle chain. The elongate flexible member may be a 
belt made of materials including, rubber, plastic, and leather. The elongate flexible member may have a plurality of 
metal elements spaced throughout the belt. 

[0025] Adverting to all of the drawings, the apparatus of the invention broadly comprises vibration sensor 24, speed 
sensor 22 and computer module 35. The vibration sensor measures the vibration of chain 18; the speed sensor meas- 
45 ures the speed of the chain and the computer module calculates power output, which power output is calculated from 
the tension and speed values. The power output value is displayed on a LCD display which is integral with the computer 
module. 

[0026] Chain speed is measured as the chain passes through the rear derailleur. As shown in Figure 2, derailleur 
19 comprises cage 50, pulley wheel 30 and jockey wheel 29. Both the pulley wheel and the jockey wheel are arranged 

50 for rotation within the cage. Pulley wheel 30 rotates about mounting bolt 51, whereas jockey wheel 29 rotates about 
mounting bolt 36. Chain 18 engages both pulley wheels and is operatively arranged to rotate pulley wheel 30 in a 
counterclockwise direction and jockey wheel 29 in a clockwise direction when the chain is driven in a clockwise rotation 
by front sprocket 16. Jockey wheel 29 comprises a plurality of embedded magnets 29a, 29b, 29c, 29d, and 29e. The 
magnets are AINiCo in composition and are equally spaced radially and angularly. Although five magnets are shown 

55 in the drawing, any number of magnets can be used. In fact, the attached software code for the invention is written for 
four equally spaced magnets. 

[0027] Speed sensor 22 comprises approximately 4000 turns of 40 gauge enamel coated wire, wound about a 0.1 25 
inch diameter mild carbon steel dowel, which is then bent to form a U-shape. The winding has a resistance of approx- 
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imately 200 ohms. Once the winding is made, it is potted in place with epoxy. As shown in Figures 2 and 3, sensor 22 
is mounted to bracket 50 by bolt 36. The body of the sensor may have an integral mounting bracket 38 or may be 
mounted using any standard mounting means known in the art. The sensor may be a proximity sensor such as a simple 
reed switch or some other type of sensor including, but not limited to, magnetic, inductive, or capacitive sensors. 
s [0028] In operation, speed sensor 22 detects the movement of the magnets as pulley wheel 29 rotates. For a known 
pulley diameter, the linear speed of the chain can be derived from the measured rotational speed of the pulley wheel 
from the following equation: 



10 



v = cor 



where 

v = linear chain speed; 
75 co = rotational chain speed; and 

r = radius of the pulley wheel 

[0029] Although one means of measuring chain speed has been described in detail, it should be readily appreciated 
by those having ordinary skill in the art that other means of measuring chain speed are possible. For example, a magnet 

20 (or magnet ring containing multiple magnetic poles) can be mounted to one of the bicycle cranks. A speed sensor can 
be mounted to the bicycle's wheel. Crank speed, wheel speed, and a table of the bicycle's gear ratios then allows the 
calculation of chain speed. Alteratively, a magnet can be mounted to one of the cranks and the sensor can be a gear 
position indicator. This indicator would sense which of the bicycle's front sprockets is engaged with the chain. These 
sensors, along with a table of the bicycle's gear ratios then allows the calculation of chain speed. Finally, chain speed 

25 can be detected directly using a magnetic sensor, such as a variable reluctance (VR) sensor, a Hall Effect sensor or 
a magneto resistive sensor. This method measures the linear movement of chain pins past the sensor to produce a 
signal that is representative of chain speed. 

[0030] The next step in determining power is to measure the vibration or tension of the drive chain. Examples of 
sensors that can be used to measure the vibration include but are not limited to variable reluctance sensors, Hall effect 

30 sensors, magneto resistive sensors, inductive, capacitive, optical and acoustic sensors. These sensors measure the 
frequency of vibration of the drive chain. In the case of a ferrous object, such as a bicycle chain moving past a magnetic- 
type sensor such as a variable sensor, an electric current is generated that is characteristic of the body's vibration. 
Such vibration results from excitation, or "plucking", which is a by-product of the normal operation of most devices, 
due to the inevitable friction and slight roughness of the tensioned body engaging with its drive or guide interfaces, 

35 such as rollers, sprockets, etc. In the case of a bicycle, this vibration is caused by the constant lateral and longitudinal 
accelerations of the bicycle, the motion induced by the individual rollers of the chain engaging and disengaging with 
the teeth of the drive sprockets, and the roughness of the surface on which the bicycle is being pedaled. 
[0031] In a multi-speed bicycle the chain is capable of being shifted to a number of different sprockets at both the 
front and the rear. This causes the chain's position (relative to the bike) to change, but the overall length of the vibrating 

40 section remains essentially constant. This is because the chain engages drive sprockets at both ends of the tensioned 
section, and the position of the axis of rotation of these sprockets does not change. The constant length likewise means 
the mass of the vibrating section remains essentially constant. 

[0032] In a preferred embodiment, vibration sensor 24 is a variable reluctance sensor (VR sensor). As shown in 
Figures 4 and 5, the VR sensor comprises a 1 .0" x 0.5" x 0.25" rectangular A1 NiCO magnet 40 wrapped with approx- 

45 imately 3000 turns of 40 gauge copper enamel coated wire having a resistance of 1 .049 ohms/foot. The total resistance 
of the coil is approximately 1 000 ohms. The faces of the magnet are laminated with rectangular carbonite plastic sheets. 
The winding is potted into place with epoxy. The sensor is housed in a plastic casing 42 and is mounted to chain stay 
23, such that its face 44 is in close proximity to the upper run of chain 18. In a preferred embodiment, as shown more 
clearly in Figures 5 and 6, the sensor is mounted to the chain stay by means commonly known to those skilled in the 

50 art. so that face 44 of the sensor is directly below the upper run of the chain. During the operation of the bicycle, as 
the chain is moved from sprocket to sprocket in the front and rear, thereby changing the gear ratio, the position of the 
chain changes depending upon the front and rear sprockets it is on. As shown in Figure 6, the sensor is mounted at 
an angle so that face 44 of the sensor is at a constant distance from the chain irrespective of the sprockets being used 
by the rider. In a preferred embodiment, this angle is approximately 30 degrees. Nevertheless, a change in the distance 

55 between the sensor and the chain (due to shifting gears, for example) will affect the amplitude of the signal, but not 
the frequency. Since the tension in the chain is related to the frequency of vibration of the chain, and not the amplitude 
of vibration, the magnitude of excitation is not relevant to the tension measurement. 

[0033] The natural frequency of vibration of a body such as a string or chain is related to the tension carried by the 
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body. This relationship is defined by the equation: 

F = 4mlv 2 

where: 

F = force carried by body in Newtons 
m = mass of tens ion ed section of body in kilograms 
I = length of tensioned section of body in meters 
v = frequency of vibration in hertz. 

[0034] In addition to the primary vibration that characterizes the tension in a body, other structures within the bicycle 
chain may cause extraneous signals. One such extraneous signal is caused by the pins of the bicycle drive chain 
moving past the. sensor. Another extraneous signal is characteristic of the chain disengaging from the sprocket teeth. 
A yet another extraneous signal is caused by the uneven distribution of the mass of the chain between the links and 
pins. In addition, other extraneous signals may be generated due to the relationship of the sensor to the chain. These 
extraneous signals can be filtered out during signal processing so as to leave only the primary vibration signal. These 
signals can be eliminated by various means. One method is to construct the sensor so that it combines the signals 
from two or more physical locations in such a way that the extraneous signal components cancel, or to use a sensor 
that has a relatively low response in the frequency band where the extraneous signal components occur. A second 
method is to use analog circuitry in the computer to filter out the frequency range of the extraneous components, which 
can be a fixed range, or a variable range controlled by the chain speed sensor. In the preferred embodiment, a third 
method is employed, in which the extraneous signal components are filtered out by means of digital signal processing, 
as described in the attached computer code for the DSP 

[0035] As previously described, the overall length of the vibrating section of the chain remains essentially constant 
because the axes of rotation of the sprockets at both ends of the chain section do not change. However, the length of 
the section of the chain that vibrates does not remain exactly constant as various front and rear sprockets are selected. 
If "the sprockets being used are differing sizes, the length of the vibrating length wilt be reduced. Depending* on the 
level of accuracy desired in the power measuring system, this variation can be ignored, as in many circumstances, the 
error introduced will be insignificant. In those cases, it is sufficient to perform calculations using a constant value for 
the length of the vibrating section of the chain. This value can be either the distance between the centers of the sprock- 
ets, or that value reduced by a small amount. 

[0036] If greater accuracy is desired, or if there will be large variation in the sizes of the sprockets used, the exact 
vibrating length can be calculated. The vibrating length is, in mathematical terms, the length of a line segment that is 
externally tangent to the circles formed by the two sprockets. In mathematical terms, the length of such a segment is: 

L = D * sqrt (1 - ((R1 - R2) / D)**2) 

where: 

• L = the desired length of the vibrating section 
D = the distance between the centers of the two sprockets (this value does not change as the gears are shifted) 
Rt = the radius of one of the sprockets 
R2 = the radius of the other sprocket 

[0037] The radius of a sprocket can be calculated from the number of teeth on the sprocket by means of the formula: 

R = T*P/(2*PI) 

where: 

R is the desired radius 
T is the number of teeth 

P is the pitch (distance between adjacent teeth) 
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[0038] The sizes of the sprockets can be determined in various ways. If the shifting mechanism provides information 
as to what sprocket from a set has been selected, or if sensors are fitted to determine which sprockets the chain is 
actually on, then the sprocket sizes can be determined from a list of available sprockets. 

[0039] Alternatively, if the chain linear speed, and the rotational speed of the two sprockets are known, sprocket 
sizes can be determined without prior information about what sprockets are available. The rotational speeds of the 
sprockets are commonly measured on bicycles, with the speed of the front sprocket used to calculate cadence, and 
the speed of the rear sprocket (which is equivalent to the speed of the bicycle wheel) used to calculate bicycle speed. 
This determination of the sprocket is by means of integrating the chain speed over one rotation of the desired sprocket. 
In a preferred embodiment, the chain speed sensor will detect links of the chain directly, and the links can be counted 
over the period of one revolution of the sprocket in question. To reduce the errors due to rounding the links can be 
counted over two or more revolutions and divided by the number of revolutions. Alternatively, if the size of one sprocket 
has been determined, the size of the other sprocket may be obtained by multiplying the size of the first sprocket by 
the ratio in speeds between the two sprockets. 

[0040] The vibration sensor and the speed sensor are electrically connected to the computer module, which contains 
a signal processing unit. The module contains a liquid crystal display (LCD) and keyboard. The module is typically 
mounted on the handlebars for easy viewing during cycling, although it can be mounted in any convenient location. In 
a preferred embodiment, the computer module comprises a Texas Instruments TMS320C3X Starter Kit, and an addi- 
tional circuit board which comprises the peripheral interface circuit. The basic components of the TMS320C3X DSP 
Starter Kit are the TMS320C31 DSP, a TLC32040 analog interface circuit (AIC), expansion connectors, system clock, 
parallel printer port interface, and tri-color LED. 

[0041] The operation of the electronic circuitry of the invention is best understood with reference to Figure 7, which 
illustrates a general schematic and block diagram of the invention. In this drawing, block 100 represents vibration 
sensor 24. As the chain vibrates in the vicinity of this sensor, a signal with a voltage in the range of approximately 2 
mV to 200 mV is produced. The relevant signal due to the frequency characteristic of the tension in the chain will be 
one component of this signal, but it may be accompanied by other signal components. The output signal from this 
sensor is fed to block 103, which represents the chain vibration amplifier 

[0042] Amplifier 103 amplifies the chain vibration signal to a level that is appropriate for circuit components 104 and 
105. The gain level can be adjusted two different ways, for four different gain settings. One adjustment is via jumper 
JP8, and the other is by means of the low gain signal, which is under control of the microprocessor. Although the 
amplitude of the signal from sensor 1 00 can vary, the dynamic range of circuit components 104 and 105 is high enough 
that the gain level should not need to be adjusted during operation; these two adjustments can be made once and left 
alone. The output of this section is an amplified version of its input, and is fed to circuit components 104 and 105. 
Alternatively, the output signal from this sensor is fed to block 104, which represents the low-pass filter. 
[0043] Circuit component 104 is a low-pass filter and component 105 is a 14 bit A/D converter. These two circuits 
are both contained in the Texas Instruments TLC32040 AIC part on the Starter Kit board. This part can be configured 
to convert its analog input to a 14-bit digital stream, sampled at 450 Hz, and it includes a Nyquist frequency low-pass 
filter to eliminate aliasing. The digital output of this section is read by the DSP (circuit component 110) for frequency 
calculations. 

[0044] Component 101 contains chain speed sensor 22 as described previously. The sensor comprises a wire coil 
mounted on rear derailleur cage 50. Jockey wheel 29 of the derailleur contains embedded magnets 29a, 29b, 29c, 29d 
and 29e that induce a voltage in the coil as they pass by. This results in an oscillating signal whose frequency is directly 
proportional to the rotational speed of the jockey wheel, and consequently, to the linear speed of the chain that is driving 
the jockey pulley. The amplitude of this signal is typically in the range of 20 mV to 2V. This signal is fed to hysteresis 
circuit 102. It should be noted that, although five magnets are embedded in the jockey wheel in the preferred embod- 
iment, any number of magnets can be embedded. Also, other means of measuring chain speed may be employed. 
[0045] Hysteresis circuit 1 02 converts the output signal of chain speed sensor 101 from its low voltage level to a TTL 
level that can be read by digital signal processor (DSP) 110. Intermediate level information is not of interest-onlythe 
fundamental frequency due to the passing of each jockey wheel magnet is relevant. For this reason, the circuit consists 
of a Schmitt Trigger, which has very high gain to "clip" the signal, and which has hysteresis to suppress any low-level 
noise that may be present in the signal. The result is a pulse train with each pulse corresponding to one of the jockey 
wheel magnets. This signal is input to digital signal processor 110. 

[0046] Digital signal processor (DSP) 110 and random access memory (RAM) 106 operate together to perform the 
power calculations. In a preferred embodiment, the DSP is a Texas Instruments Model No. TMS320C31. The RAM 
contains the program code for this device (although such code could also be stored in a ROM chip connected to the 
DSP in an analogous manner). The DSP calculates power based on two inputs. The first is the signal from circuit 
component 102, which can be used to determine the linear speed of the chain (given the fixed values of how many 
magnets are on the jockey wheel, and how many teeth are on the jockey wheel). The second is the signal from circuit 
components 1 04 and 1 05, which can be used to determine the vibrational frequency of the chain, and in turn the tension 
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of the chain (the speed of the chain, determined from the jockey wheel signal, is also used to assist in selecting the 
proper component from the output signal of circuit components 104 and 105). The tension of the chain is multiplied by 
the speed of the chain to give a result of the power transmitted through the chain at any instant. The method used to 
do this calculation is described fully in the attached computer code for the DSR 

[0047] The signal is processed using a standard technique known to those skilled in the art including, but not limited 
to, Fast Fourier Transform, Discrete Fourier Transform, Fourier Transform, Zero-Crossing, 1 -bit FFT with cosine/sine 
table look-up, wavelets, filter banks, FIR filters, MR filters, adaptive filter techniques (including Least Mean Square 
Algorithm), adaptive interference cancelling, adaptive self-tuning filters, adaptive line enhancers, and gabor expansion. 
In one embodiment, particularly if there are no significant interference signals present, the frequency of the vibration 
can be determined by using "zero-crossing" algorithm, which measures the period of a sinusoid wave by realizing when 
the signal passes through the median voltage and measuring the time period between median voltage crossings. In 
another embodiment, a Fast Fourier Transform (FFT) is performed on the signal. Linear prediction algorithm (LPA) and 
auto regression method (AR) can be used to find more practical and more implementable solution in a microprocessor 
based technique. 

[0048] To decrease the amount of calculation capacity for the vibration frequency determination the chain speed 
variation signal can also be used. Because the chain speed reflects the power and force variation against the pedals 
and cranks this chain speed variation signal can be used as a reference signal form for the power variation signal form. 
And because these signals are closely correlating it is enough to calculate absolute power only sometimes so that the 
absolute power level of varying signal can be fixed. So it is necessary calculate and define the absolute power of force 
and tension only some (for example 1-2 ) times per one cadence cycle. That is because the power signal form is now 
known from the chain speed signal and some points are used to fix that signal to the absolute scale values. 
[0049] In a preferred embodiment of the invention the speed of the elongate flexible member can be measured 
directly by using a speed sensor . In another preferred embodiment of the invention the speed is determined indirectly, 
for example by using cadence measurement. In said another embodiment the speed of the elongate flexible member 
is determined by using cadence signal, wheel speed and sprocket sizes or the ratio of sprocket sizes. 
[0050] Chain speed is calculated by counting the time between the rising edges of the jockey wheel signal. 
[0051] Memory 107 is a non-volatile memory chip, such as a serial EEROM. The DSP is operatively arranged to 
write and read data to and from memory 107. Memory 107 is used to store installation- specific information such as 
chain-line length and chain linear density. 

[0052] Circuit components 108 (liquid crystal display (LCD)), 109 (keyboard/CPU), and 111 (keyboard) comprise a 
simple "terminal" operating at TTL voltage levels and typically arranged to be mounted on the handlebars of the bicycle. 
The DSP can send the calculated power value (or any other relevant information) to this CPU, which will display the 
information on the LCD. Key presses on the keyboard will result in information being sent to the DSP, which can respond 
with some appropriate action (e.g., changing the displayed power value from watts to horsepower.) Any suitable user 
interface can be implemented. These three components comprise a standard off-the-shelf unit manufactured and avail- 
able by Bolton Engineering of Melrose, MA. In a preferred embodiment, a Bolton unit which uses a Microchip PIC16C63 
processor was used. Of course, other processor units would also function satisfactorily In a preferred embodiment of 
the invention the computer generates a signal representing the power output of the elongate flexible member to a 
external unit coupled to said computer, said external unit having a display for displaying the power output. Said external 
unit is the terminal referred above. In a preferred embodiment of the invention the external unit coupled to said computer 
is a heart rate monitor having a display. Heart rate sensor measures the ECG of the bicycle user. 
[0053] Circuit component 112 is an optional RS232 interface, well known to those having ordinary skill in the art. 
This interface converts the TTL signal levels of the UART connected to the DSP into RS232 levels for communication 
with an external device such as a personal computer (PC). The DSP can use this interface to upload a series of saved 
power values to a host PC for analysis. 

[0054] Component 113 is a host PC. This PC is an external device which is not used in normal operation of the 
invention. It can be any sort of computer that can retrieve information from the power meter device for further processing. 
It would typically be connected to the power meter at the completion of a bicycle fide. 

[0055] Figures 8-11 show electrical schematic diagrams of the peripheral interface circuit of the invention. The pe- 
ripheral interface circuit connects the several peripheral devices (e.g., RAM, keyboard, display, etc.) of the invention 
to the microprocessor unit 110. The signal connections in the physical connectors of the electronic unit are shown in 
Fig. 8 and the connections of a keyboard and display in Fig. 9 and the connection of a memory unit in Fig. 10. The 
connection of the sensors such as wheel speed sensor (WHEEL), cadence sensor (CADENCE), chain speed from the 
jockey wheel (JOCKEY) and vibration sensor (AUX) are shown in Fig. 11 to the Texas Instruments TMS320C3X DSP 
Starter Kit board. 

[0056] Figure 1 2 is a flowchart of the computer code used to run the DSP in a preferred embodiment of the invention. 
The code itself is as follows: 
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dsp.c - processing code for powermeter monitoring system 



Copyright 1997 - CC Kinetics, All Rights Reserved 



[0057] This code must be called every 2 seconds with BLOCK_SIZE number of struct sJockAD(s). Each struct has 
a member for a/d data and jockey zero-crossing data. The a/d data must be in integer offset notation, where the offset 
is specified by the macro AD_OFFSET. The jockey data is a one-bit value, where a crossing is represented by a 1 and 
*5 a non -crossing by zeros. 

[0058] Invoke the function ProcessPowerlnit to initialize the processing code. Call the function ProcessPower every 
2 seconds with a/d and zero crossing values as the first argument. Results are returned through the function's second 
argument. 

20 

#include<stdio.h> 
25 #include<math.h> 

#include<malloc.h> 
#include<string.h> 



30 



#define RAW_SAMP_FREQ 7200 
#define RAW_SAMP_PERIOD 
((float) 1 /(float)RAW_SAMP_FREQ) 
35 #define RESAMPLE_FACTOR 1 6 

#define FFT_SIZE 64 
#define FFT_OFFSET 4 
40 #define BLOCK_PERIOD_SEC 2 

#define BLOCK_SIZE 
(RAW_SAMP_FREQ * 
BLOCK_PERIOD_SEC) 
#define BLOCK_SIZE_RESAMPLE 
(BLOCK_SIZE/RESAMPLE_FACTOR) 

so 
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#define N U M_FFT_PE R_B LOC K ((BLOCK_SIZE_RESAMPLE- 
FFT_SIZE+ 

FFT_OFFSET)/FFT_OFFSET) 

#define AD_OFFSET 0x80 
#define NUMJNTERVALS 3 
#define JOCK_REV_RATIO 4 
#define JOCK_NUM_TEETH 10 
#define PINS PER METER 



((float)78.7) 



#deflne M_PI 

3.14159265358979323846 
#define M_2PI 
(2.0*M_PI) 
#define SQ(x) 
((x)*(x)) 

#define MIN_RPM 
((float)20.0) 
#define MAX RPM 
((float)200.0) 

#define MIN_CHAINRING 

#define MAX_CHAINRING 

#define MIN_PIN_FREQ 
MIN_CHAI N RING/(float)60) 

#define MAX_PIN_FREQ 
MAX_C HAINRIN G/(float)60) 

#define MIN_JOCK_FREQ 
* JOCK_REV_RATIO/ 

JOCK_NUM_TEETH) 

#define M AX_J OC K_F RE Q 

(MAX_PIN_FREQ * JOCK_REV_RATIO/ 

JOCK_NUM_TEETH) 

#define PAD_SIZE 
((int)(RAW_SAMP_FREQ/MIN_JOCK_FREQ)) 



((float)24.0) 
((float)60.0) 



(MIN_RPM 
(MAX_RPM 



(MIN_PIN_FREQ 



#define BINS_PER ((FFT_SIZE/2)+1) 
#define SAMP_FREQ 
(RAW_SAM P_F REQ/ 
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(float)RESAMPLE_FACTOR) 
#define SAMP_PERIOD 
((float) 1 /(float)SAMP_FREQ) 
#define NUM_FFT_PER_SEC 
(SAMP_FREQ/FFT_OFFSET) 
#define WINDOW_L0BE_WIDTH 
#define INTERVAL_1_SEC 
#define INTERVAL_2_SEC 
#define INTERVAL_3_SEC 
#define NUM_FFT_PER_1_AVE 
NUM_FFT_PER_SEC)) 
#define N U M_FFT_P ER_2_AVE 
NUM_FFT_PER_SEC)) 
#define NUM_FFT_PER_3_AVE 
NUM_FFT_PER_SEC)) 
#define MAX_AVE RAGE_WI N DOW 
((inf)NUM_FFT_PER_3_AVE) 
#define MAX_FILTER_COEF 
#define BIN2FREQ(b) 
(((float)b)/((float)FFT_SIZE/ 
SAMP_FREQ)) 
#define FREQ2BIN(f) 
((float)FFT_SIZE/SAMP_FREQ)) 

#define PINHZ2SPEED(f) 
#define ATTEN_NOTCH_DB 
((float)30) 
#define MM2M 
((float) 1000) 
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((int)(INTERVAL_1_SEC * 
((int)(INTERVAL_2_SEC * 
((int)(INTERVAL_3_SEC * 



41 



((f) 



((f)/PINS_PER_METER) 



extern float fChainLength_mm; 
extern float fChainLinearDensity; 
extern Int nChainring; 

struct sPowerAve { 

int nTwang[NUM_INTERVALS]; 
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int n P in Hz[N U MJ NTERVALS]; 
int nPower(NUMJNTERVALS]; 



}; 



10 struct sPowerArray { 

float fTwa ng [N U M_F FT JP E R_B LOC K] ; 
float fPi n Hz[N U M_F FT_P E R_B LOCK]; 
float fPower[NUM_FFT_PER_BLOCK]; 



}; 



struct sPowerSet { 
struct sPowerAve ave; 
struct sPowerArray array; 

}; 

struck sJockAD { 
unsigned int jock: 1; 
unsigned int ad:31; 

}; 

struct sCmplex { 
float re.im; 

}; 



45 struct sFiltCoef{ 

int num; 



float coef[MAX_FILTER_COEF]; 

}; 

struct sFiltCoefHistory { 
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struct sFiltCoef flit; 

float histoiy[MAX_FILTER_COEF]; 

}; 

struct sAverageFilter { 
int length; 
int index; 
int count; 

float value[MAX_AVERAGEJ/VINDOW]; 

}; 

static struct sComplex cData[FFT_SIZE]; 

static struct sAverageFilter 

AveTwang[NUMJNTERVALS],AvePinHz 
[NUMJNTERVALS],AvePower[NUMJNTERVALS]; 

static struct sFiltCoefHistory decFilter; 

static int nFftCount; 

static int nBlockCount; 



0909940A2 I > 



13 



EP 0 909 940 A2 



»********************************»***««***«»*»***»******»».*»»** 



FilterLP4000-FIR filter, LP=300hz, direct realization, 30db atten. 

static struct sFiltCoefFilterLP4000={ 
41, 

-0.005125881425025383, 

0. 002582302 1 50794436 , 

0.01 052733850 1 65874. 

0.01681942643599849, 

0.01973515196700159. 

0.01812679354307367, 

0.01 174820956906760, 

0.001421601768124991, 
-0.01 1 00256787642489, 
-0.02290073789531828, 
-0.03131135456545855, 
-0.03347630998683525, 
-0.02738957008440296. 
-0.01224700060318999, 

0.01 129467944240988. 

0.04110924668809402. 

0.07385658812070439. 

0.1054668790709014, 

0.1317734248780911, 

0.1491843151328381, 

0.1552734374999997. 

0.1491843151328384, 

0.1317734248780913. 

0.1054668790709016, 

0.07385658812070449. 

0.04110924668809407, 

0.01129467944240987, 
-0.01224700060319003, 
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-0.02738957008440303, 
-0.03347630998683532, 
-0.03131135456545862, 
-0.02290073789531 832, 
-0.01100256787642492, 
0.001421601768124995, 
0.01174820956906763, 
0.01812679354307370, 
0.01973515196700161, 
0.01681942643599851, 
0.01052733850165875, 
0.0025823021 50794424, 
-0.005125881425025422 

}; 

y*** *********** 

log2 

*********** 

static int log2(unsigned int x) 
{ 

unsigned int mask, I; 
if (x=0) 

return- 1; 
x~; 

mask=1; 

l=0; 
while(1){ 
if (x=0) 

retum(i); 

x=x&(-mask); 

mask *=2; 
i++; 
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windowGaussian 

**********^***************************************^***#******* 

static void windowGaussian(float dataQ, int N, int a) 
{ 

int n; 

static int Nstore=0; 
static int Astore=0; 
static float loAlpha; 
static float *W; 
if (Nstore!=NIIAstorei=a){ 

if (Nstore) 

free(W); 

Nstore=N; 

Astore=a; 

W=(float *)malloc(sizeof(float) * N); 

for (n=0; n<N; n++) 

data[n] *=W[n]; 

} 



fft 

*^*************«**********^**#******«*********«**************** 

static void fft(struct sComplex *x, int mm) 
{ 

static struct sComplex *w; 
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static int mstore=0; 
static int n=1; 

struct sComplex u,temp,tm; 
struct sComplex *xi,*xip t *xj,*wptr; 

int ij.k.Ue.windex; 

int m=log2((unsigned int)mrn); 

float arg , w_re , vvjm.wrecurjm ,wtemp_re ; 

if(m!=mstore){ 
if(mstore)!=0) 
free(w); 
mstore=m; 
if(m==0) 

return; 

n=1«m; 
le=n/2; 

w=(struct sComplex *)calloc(le-1 ,sizeof(struct sComplex)); 
if (!w) 

return; 

arg=(float) (M_PI/I); 
wrecur_re=w_re=float)cos(arg); 
wrecur_im=wJm=float)-sin(arg); 
xj=w; 

for 0=1;j<le;j++H 

xj->re=(float)wrecur_re; 
xj-> i m= (float) wrecu r j m ; 
xj++; 

wtemp_re=wrecu r_re*w_re- 

wrecurjm*w_im; 

wrecurjm=wrecur_re*wjm+wrecurjm*w_re; 

wrecur_re=wtemp_re; 
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} 

} 



le=n; 

windex=1; 
for(1=0;1<m;1++){ 
le=le/2; 



for(i=0;i<n;i=i+2*le){ 
xi=x+i; 
xip=xi+le; 

temp.re=xi->re+xip->re; 
temp.im=xi->im+xip->im; 
xip->re=xi->re-xip->re; 
xip->im=xi-> im-xip-> im ; 
*xi=temp; 



wptr=w+windex-1 ; 
forG=1;j<le;j++){ 

u=*wptr; 

for(i=j;i<n;i=i+2*le){ 
xi=x+i; 
xip=xi+le; 

temp.re=xi->re+xip->re; 

tempjm=xi->im+xip->im; 

tm.re=xi->re-xip->re; 

tmJm=xi->im-xip->im; 

xip->re=tm.re*u.re4m.im*u.im; 

xip->im=tm.re*u.im+tm.im*u.re; 

*xi=temp; 

} 

wptr=wptr+windex; 

} 

windex=2*windex; 

} 
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j=0; 

for(i=1;i<(n-1);i++){ 
k=n/2; 
while(k<=j){ 

j=j-k: 
k=k/2; 

} 

j=j+k; 
if(i<j){ 

xi=x+i; 

xj=x+j; 

temp=*xj; 

*xj=*xi; 

*xi=temp; 



} 



} 
} 



AverageWindowlnit 



static void AverageWindowlnit(struct sAverageFilter *af, int 
40 numPoints) 

{ 

if(numPoints>MAX_AVERAGE_WINDOW) 

exit(-1) 
af->count=af->index=0; 
af->length=numPoints; 

} 



55 



19 



NSOOCID: <EP 0909940A2 I > 



EP 0 909 940 A2 



AverageWindowUpdate 
**********************************^*** 

static float AverageWindowUpdate(struct sAverageFilter*af, float 

value) 

{ 

int i; 

float ave; 

/*put value in filter 
V 

af->value[af->index]=value; 
af->index=(af->index+1)%af->length; 
if(af->count<af->length) 
++af->count; 

^average it 
*/ 

ave=0; 

for(i=0;i<af->count;i++) 
ave+=af->value[i]; 

return ave/af->count; 
} 

FIRfilterlnit 

static void FIRfilterlnit(struct sFiltCoefHistory *fir) 
{ 

int i; 
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for(i=0;i<fir->filt.num;i++) 

fir->history[il=0.0; 

} 



FIRfilterSamp 

************** ****** ********* *********************************** j 

static float FIRfilterSamp(float input, struct SFiltCoefHistory*fir) 
{ 

int i; 

float *histj3tr,*hist1_ptr,*coef_ptr; 
float output 

hist_ptr=fir->history; 
hist1_ptr=hist_j)tr; 
coef_ptr=fir->filt.coef+fir->filt.num-1 ; 

output=*hist_ptr++*(*coef_ptr~); 
for(i=2;i<fir->filt.num;i++){ 

*hist1_ptr++=*hist _ptr; 

output+=(*hist__ptr++)*(*coef_ptr~); 

} 

output+=input*(*coef_ptr); 
*hist1_ptr=input; 

return output; 

} 

I**************************************************************** 
FIRdecirnate 

****************************************************************/ 
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static void FIRdecimate(float *input, float *output, int N, struct 
sFiltCoefHistory*fir, int decimateFactor) 

{ 

int ij; 

for(i=0;i<N;i++) 
input[i]=FIRfilterSamp(input[i] l fir); 

for(i=O t j=0;i<N;i+=decimateFactorj++) 
output[j]=input[i] 

} 

GetADvalues 

static void GetADvalues(struct sJockAD *pJockAD, float *pAD) 
{ 

int i; 

for(i=0;i<BLOCK_SIZE;i++) 
pAD[i}=((float)pJockAD[i].ad-AD_OFFSET); 

} 

^************ ******************* ****************^ 
Decimate 

****************************************************************^ 

static void Decimate(float *indata, float *outdata) 
{ 
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FIRdecimate(indata, outdata, BLOCK_SIZE, &decFilter, 
RESAMPLE_FACTOR); 

} 

ComputePinHz 

**************************** ******************************** 

static void ComputePinHz(struct sJockAD *pJockAD, float *pPinHz) 
{ 

float xcenter, x1 , x2, y1 , y2 t period; 
float xlefM. xleft2, xrightl, xright2, m, b; 
int i,j, found; 

for (i=0;i<NUM_FFT_PER_BLOCK;i++){ 

xcenter+(i * RESAMPLE_FACTOR * FFT_OFFSET)+(FFT_SIZE * 
RESAMPLE_FACTOR/2)+PAD_SIZE; 

Xleft1=xleft2=xright1=xright2=0; 

found=0; 

for(j=xcenter;j>0;j~) 

if(pJockAD[jljock){ 

found=1; 
break; 

} 

if(found) 

xleft1=j; 

else{ 

pPinHz[i]=0; 
continue; 
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} 



found=0; 

forG=xcenter+1;j<BLOCK_SIZE+2*PAD_SIZE;J++) 
if(pJockAD[j].jock){ 
found=1; 
break; 

} 



if(found) 
else{ 



xright1=j; 

pPinHz[i]=0; 
continue; 



} 



if(xcenter-xleft1 <xright1 -xcenter){ 
found=0; 

forQ=xlefl1-1;j>0;H 

if(pJockAD[j].jock){ 
found=1; 
break 

} 



if (found) 

xleft2=j; 

else{ 

pPinHz[i]=0; 
continue; 



} 

x1=(xleft1+xleft2)/2; 
x2=(xright1+xleft1)/2; 
y1=xleft1-xleft2; 
y2=xright1-xIeftT; 
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} 

else{ 



} 



found=0; 

for(j=xright1 +1 ;j<BLOCK_SIZE+2*PAD_SIZE;j++) 
if(pJockAD[j].jock){ 
found=1; 
break; 

} 

if(found) 

xright2=j; 

else{ 

pPinHz[i]=0; 
continue; 

} 

x1=(xright1+xleft1)/2; 
x2=(xright1+xright2)/2; 

y1-xright1-xleft1; 
y2-xright2-xright1; 



m=(y2-y1)/(x2-x1); 
b=y2-m * x2; 
period+m * xcenter+b; 

pPinHz[i]=((RAW_SAMP_FREQ/period)JOCK_REV_RATIO)* 
JOCK_NUM_TEETH 
} 
} 



****************************************** 
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NotchFilter 

static voidNotchFilter(float*x, float freq, int N) 
{ 

const int sigma=8; 

const int nHzStart=(int)MIN_PIN_FREQ; 

const int NHzStop=(int)iyyiAX_PIN_FREQ; 

const int nNumVectors=nHzStop-nHzStart+1; 

const float fMinNotch-0.01f; 

const float fAmpNotch-100.f; 

static int Nstore=0; 

static struct sComplex *carNotch; 

static float *fNotch[(int)(MAX_PIN_FREQ-MIN_PIN_FREQ+1)]; 
int Nfft=(n-1)*2; 
int ij; 

float f,fMax,atten; 
if (N!=NStore){ 
if(NStore|=0){ 

free(carNotch); 

for(i=nHzStart; l<nHzStop;i++) 
free(fNotch(i)); 

} 

if (N==0) 
return; 

carNotch=(struct sComplex *)calIoc(Nfft, 
sizeof(struct sComplex)); Nstore=N; 

for(i=0;i<nNumVectors;i++){ 

fNotch[i]=(float*)calloc(Nfft,sizeof(float)); 

for(j=0;j<Nfft;j++) 
fNotch[i]0]=fAmpNotch*(float)sin(M_2PI*(fIoat) 
(i+nHzStart) * SAMP_PERIOD * j); 
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windowGaussian(fNotch[i], Nfft, sigma); 

for(j=0;j<Nffi;j++){ 

carNotch[j].re=fNotch[i][j]; 
carNotch[j].im=0; 

} 

fft(carNotch, Nfft); 

forO-0;j<N;j++) 

fNotch[i]D]=SQ(carNotchD]re)+SQ 

(carNotchOl-im); 

fMax=fNotch[iK0]>fMinNotch/2?fNotch[i][0] : 

fMinNotch/2; 

for(j=1 ;j<N;j++) 

if(fNotch[i][j]>fMax) 

fMax=fNotch[i]01; 

attn=(float)pow(1 0. ATTEN_NOTCH_DB/(float)1 Oy 
SQ(FFT_SIZE); 

for(j=0;j<N;j++){ 

f=fNotch[i]Q]; 

f=-f; 

f/=fMax; 

f*=(float) 1 -((float) 1 /atten); 
f+=1; 

fNotch[i][j]=f; 

} 

} 

} 

if(freq>=nHzStart&&freq<=nHzStop) 
for(j=0;j<N;j++) 
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x[fl*=fNotch[(int)freq-nHzStart][j]; 

} 

Detection 

************^*************************************#***#*********^ 

static float Detection(float *data) 
{ 

int n, bin; 
float max; 

bin=0; 
max=0; 

for(n=1; n<BINS_PER;n++) 
if(data[n]>max){ 

max=data[n]; 
bin=n; 

} 

return BIN2FREQ(bin); 
} 

Power 

static float Power(float fTwang, float fChainspeed) 
{ 

float fTension; 
float fWatts; 

ffension=(float)(4.0 * SQ(fChainLength_mm/MM2M) 



* 
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* fChainLinearDensity * SQ(fTwang)); 
Fwatts=fTension * fChainspeed; 

return fWatts; 
} 

y****************************************************^ 

ProcessPower 

************ *********************************************** *****y 

voidProcessPower(struct sJockAD *pJockAD, struct sPowerSet 
*psResults) 
{ 

static struct sJockAD arCirc[bLOCK_SIZE * 2]; 
static float arAD1 [BLOCK_SIZE]; 
static float a rA D2 [ B LO C K_S I ZE_RES AM P LE] ; 
static float arAD3[FFT_SIZE]; 
static float arAD4[FFT_SIZE]; 
static float arPinHz[NUM_FFT_PER_BLOCK]; 
static struct sComplex cAD[FFT_SIZE]; 

float fPinHz, fPower t fTwang; 
int ij.k; 

40 if(nBlockCount>0){ 

memmove(arCirc,&arCirc[BLOCK_SIZE] f si2eof(struct 
sJockAD)*BLOCK_SIZE); 

memcpy(&arCirc[BLOCK_SIZE], pJockAD, 
sizeof(struct sJockAD)*BLOCK_SIZE); 

} 

else{ 

so memset(arCirc,0,sizeof(struct 
sJockAD)*BLOCK_SIZE); 
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memcpy(&arCirc[BLOCK_SIZE], pJockAD, 
sizeof(struct sJockAD) * BLOCK_SIZE); 

} 

GetADvalues(&arCirc[BLOCK_SIZE-PAD_SIZE], arAD1 ); 
Decimate(arAD1 , arAD2); 

ComputePinHz(&arCirc[BLOCK_SIZE-2*PAD_SIZE], arPinHz); 
for(i=0, k=0; k< N U M_FFT_PER_B LOC K; i+=FFT_OFFSET, k++){ 

MEMCPY(arAD3, &arAD2[i], FFT_SIZE*sizeof(float)); 

windowGaussian(arAD3, FFT_SIZE, WINDOW_LOBE_WIDTH); 

for(j=0; j<FFT_SIZE;j++){ 

cADQ].re=arAD3[j]; 
cADD]im=0.0; 

} 

fft(cAD, FFT_SIZE); 

for (j=0; j<BINS_PER; j++) 

arAD4[j]=SQ(cAD[j].re)+SQ(cAD[j].im); 

fPinHz=arPinHz[k]; 

NotchFilter(arAD4, fPinHz,BINS_PER); 
fTwang=Detection(arAD4); 
fPower=Power(fTwang, P!NHZ2SPEED(fPinHz)); 

psResults->array.fTwang[k]=fT wang; 

psResults->array.fPinHz[k]=fPinHz; 

psResults->array.fPower[k]=fPower; 

psResults->ave.nTwang[0]=(int)AverageWindowUpdate 
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(&AveTwang[0],fTwang); 

psResults->ave.nTwang[1]=(int)AverageWindowUpdate 
(&AveTwang[1 ],fTwang); 

psResults->ave.nTwang[2]=(int)AverageWindowUpdate 
(&AveTwang[2] t fTwang); 

psResults->ave.nPinHz[0]=(int)AverageWindowUpdate 
(&AvePinHz[0],fPinHz); 

psResults->ave.nPinHz[1]=(int)AverageWindowUpdate 
(&AvePinHz[1],fPinHz); 

psResults->ave.nPinHz[2]=(int)AverageWindowUpdate 
(&AvePinHz[2],fPinHz); 

psResults->ave,nPower[0]=(int)AverageWindowUpdate 
(&Ave Powe r[0] .fPower) ; 

psResults->ave.nPower[1]=(int)AverageWindowUpdate 
(&AvePower[1],fPower); 

psResults->ave.nPower[2]=(int)AverageWindowUpdate 
(&AvePower[2] .fPower) ; 

++nFftCount; 

} 

++nBlockCount; 
} 



#define SAMPLEJHZ 400 
#define SIZE_BLOCK 400 

float ComputeFreq(int data[SIZE_BLOCK]) 

{ 

int i, j; 

static float freq = 0; 

j=-i; 

for (i=0; i<SIZE_BLOCK-1; i++) { 
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if (data[i]<0 && data[i+1 ]>=0) 
if q = -1) 

j=i; 

else { 

freq = (float)SAMPLE_HZ / (float)0 - i); 
break; 

} 

} 

return freq; 
} 



so 



ProcessPowerlnit 
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voidProcessPowerlnit(void) 
{ 

nFftCount-nBlockCount=0; 
decFilter.filt=FjlterLP4000; 
FIRfilterlnit(&decFilter); 

AverageWindowlnit(&AveTwang[0],NUM_FFT_PER_1_AVE) 
AverageWindowlnit(&AveTwang[1],NUM_FFT_PER_2_AVE) 
AverageWindowlnit(&AveTwang[2],NUM_FFT_PER_3_AVE) 
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AverageWindowlnit(&AvePinHz[0],NUM_FFT_PER_1_AVE); 
AverageWindowlnit(&AvePinHz{1],NUM_FFT_PER_2_AVE); 
AverageWindowlnit(&AvePinHz(2],NUM_FFT_PER_3_AVE); 
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AverageWindowlnit(&AvePower[0],NUM_FFT_PER_1_AVE); 
AverageWindowlnit(&AvePower[1 ], NUM_FFT_PER_2_AVE); 
AverageWindowlni^&AvePowe^j.NUM.FF^PER.S.AVE); 



} 
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Claims 

I. An apparatus for measuring the power transmitted by an elongate flexible member, comprising: 

5 a system having a first engagement means separated from a second engagement means, the elongate flexible 

member coupling the first engagement means with the second engagement means so that when the first 
engagement means is rotated by a power source the elongate flexible member moves in the direction of the 
rotation of the first engagement means, which in turn causes the moving elongate flexible member to rotate 
the second engagement means; 

10 a vibration sensor that measures the vibration frequency of the elongate flexible member; 

a sensor that directly or indirectly measures the speed of the elongate flexible member; and 

a computer that determines the power transmitted by the elongate flexible member by receiving the vibration 

measurement and the speed measurement directly or indirectly measured. 

75 2. The apparatus of claim 1 , characterized in that extraneous signals are filtered by analog signal processing using 
mechanical structure of a vibration sensor or by analog electronic circuit structures. 

3. The apparatus of claim 1 , characterized in that extraneous signals are filtered by the computer using digital signal 
processing. 

20 

4. The apparatus of claim 1 , characterized in that the apparatus has a digital display for indicating the power output 
of the elongate flexible member. 

5. The apparatus of claim 1 , characterized in that the computer generates a signal representing the power output of 
25 the elongate flexible member to a external unit coupled to said computer, said external unit having a display for 

displaying the power output. 

6. The apparatus of claim 5, characterized in that said external unit coupled to said computer is a heart rate monitor 
' having a display. 

30 

7. The apparatus of claim 1 , characterized in that the system comprises a rotateable wheel being connected to said 
the second engagement means, said elongate flexible member arranged to rotate said wheel by rotating said 
second engagement means, said first and second engagement means each having at least one sprocket. 

3S 8. The apparatus of claim 7, characterized in that the speed of the elongate flexible member is determined by using 
cadence signal, wheel speed and sprocket sizes or the ratio of sprocket sizes. 

9. The apparatus of claim 1 , characterized in that the system is a bicycle. 

40 10. The apparatus of claim 1, characterized in that the elongate flexible member is a chain. 

II. The apparatus of claim 1, characterized in that the elongate flexible member is a belt having a plurality of metal 
elements equally spaced within and throughout the belt. 

4 $ 12. The apparatus of claim 1, characterized in that the first and second engagement means each comprise various 
sized sprockets so that the elongate flexible member moves from one sprocket to another sprocket on each en- 
gagement means. 

13. The apparatus of claim 12, characterized in that for determining the power of the elongate flexible member the 
so vibrating length of the elongate flexible member is calculated by using the sizes of sprockets used. 

14. The apparatus of claim 7 and 12, characterized in that the sizes of sprockets used are defined by using cadence 
signal, wheel speed and the speed of the elongate flexible member. 
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15. A method for measuring the power transmitted by an elongate flexible member in a system having a first engage- 
ment means separated from a second engagement means, the elongate flexible member coupling the first en- 
gagement means with the second engagement means so that when the first engagement means is rotated by a 
power source the elongate flexible member moves in the direction of the rotation of the first engagement means, 
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which in turn causes the moving elongate flexible member to rotate the second engagement means; the method 
comprising the steps of: 

measuring the vibration frequency of the elongate flexible member; 
5 directly or indirectly measuring the speed of the elongate flexible member; and 

receiving the vibration measurement and the directly or indirectly measured speed measurement in a compu- 
ter, 

generating a signal representing the power transmitted by the elongate flexible member. 
io 16. The method of claim 15, characterized in that extraneous signals are filtered by analog signal processing. 

17. The method of claim 15, characterized in that extraneous signals are filtered by digital signal processing. 

18. The method of claim 15, characterized in that the power output of the elongate flexible member is displayed by a 
i5 digital display of the apparatus. 

19. The method of claim 15, characterized in that the signal representing the power output of the elongate flexible 
member is provided by the computer to a external unit coupled to the computer that displays the power output. 

20 20. The method of claim 15, characterized in that the system is a bicycle. 

21. The method of claim 15, characterized in that the elongate flexible member is a chain. 

22. The method of claim 15, characterized in that the elongate flexible member is a belt having a plurality of metal 
25 material equally spaced within and throughout the belt. 

23. A method for determining the tension in a elongate flexible member on a vehicle such as a bicycle driven by the 
elongate flexible member comprising the steps of: 

30 a) identifying a vibrating section of the elongate flexible member; 

b) providing a vibration sensor on the vehicle in proximity to the vibrating section of the elongate flexible 
member; 

c) measuring the frequency of vibration of the vibrating section of the elongate flexible member by means of 
the sensor during operation of the vehicle; and 

35 d) utilizing the frequency of vibration and the known length of the vibrating section of the elongate flexible 

member to determine the tension in the elongate flexible member. 

24. The method of claim 15 or 23, characterized in that system comprises a rotateable wheel being connected to said 
the second engagement means, said elongate flexible member arranged to rotate said wheel by rotating said 

40 second engagement means, said first and second engagement means each having at least one sprocket 

25. The method of claim 1 5 or 23, characterized in that the first and second engagement means each comprise various 
sized sprockets so that the elongate flexible member moves from one sprocket to another sprocket on each en- 
gagement means. 

45 

26. The method of claim 15 or 23, characterized in that the speed of the elongate flexible member is determined by 
using cadence signal, wheel speed and sprocket sizes or the ratio of sprocket sizes. 

27. The method of claim 24 and 25, characterized in that the sizes of sprockets used are defined by using cadence 
50 signal, wheel speed and the speed of the elongate flexible member. 

28. The method of claim 23, characterized in that the length of the vibrating section of the elongate flexible member 
is identified by the known distance between the crank axle and wheel axle, said distance being input to the com- 
puter. 

55 

29. The method of claim 23, characterized in that the length of the vibrating section of the elongate flexible member 
is computed by measuring the rear sprocket size, the front gear size, and the distance between the crank axle and 
wheel axle. 
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30. The method of claim 23, further including eliminating spurious information due to links of the elongate flexible 
member moving past the sensor in the direction of the length of the elongate flexible member 

31. The method of claim 23, characterized in that said step of determining the tension in the elongate flexible member 
is performed by a computer on the bicycle. 

32. The method of claim 23, further including eliminating spurious information due to noise sources. 

33. The method of claim 23, further including eliminating spurious information due to links of the elongate flexible 
member and teeth of sprockets engaging the elongate flexible member. 

34. A method for determining the vibrational frequency of a elongate flexible member on a vehicle such as a bicycle 
driven by the elongate flexible member comprising: 

a) providing the vibration signal of the elongate flexible member by means of a vibration sensor in proximity 
to the elongate flexible member; 

b) performing a Fast Fourier Transform on the vibration signal of the elongate flexible member; and 

c) heuristically selecting the vibration frequency of the elongate flexible member 

35. The method of claim 34, characterized in that said steps of performing a Fast Fourier Transform and selecting the 
vibration frequency of the elongate flexible member are performed by a computer on the vehicle. 

36. The method of claim 34, characterized in that Linear prediction algorithm (LPA) or auto regression(AR) method is 
performed to analyse the vibration signal of the elongate flexible member. 

37. Apparatus for measuring the tension in a elongate flexible member on a vehicle such as a vehicle driven by the 
elongate flexible member comprising: 

a) a vibration sensor mounted on the vehicle in proximity to a vibrating section of the elongate flexible member 
for measuring the frequency of vibration of the vibrating section of the elongate flexible member during oper- 
ation of the vehicle; 

b) signal processing means on the vehicle operatively connected to the vibration sensor for selecting a meas- 
ured frequency of vibration; and 

c) computer means on the vehicle operatively connected to the signal processing means for computing the 
tension in the elongate flexible member 

38. The apparatus of claim 37, characterized in that the signal processing means includes means for eliminating spu- 
rious information signals due to noise. 

39. A method for measuring the power output of a person operating a vehicle such as a bicycle driven by means 
including a elongate flexible member comprising the steps of: 

a) measuring the vibration of the elongate flexible member by means of a sensor on the vehicle; 

b) measuring directly or indirectly the speed of the elongate flexible member by means of a sensor on the 
vehicle; 

c) computing the power transmitted by the elongate flexible member by means of a computer on the vehicle 
using the directly or indirectly measured speed of the elongate flexible member and the vibration of the elongate 
flexible member 

40. The method of claim 39, further including displaying the computed power by means of a display. 

41 . The method of claim 39, further including eliminating spurious information due to noise sources by means of signal 
processing means. 

42. The method of claim 41 , characterized in that said step of eliminating spurious information comprises signal 
processing employing Fast Fourier Transform. 

43. Apparatus for measuring the power output of a person operating a vehicle such as a bicycle driven by means 
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including a elongate flexible member comprising: 

a) a vibration sensor mounted on the vehicle in proximity to a vibrating section of the elongate flexible member 
for measuring the vibration frequency of the elongate flexible member; 
5 b) a sensor mounted on the vehicle in proximity to the elongate flexible member for directly or indirectly meas- 

uring the speed of the elongate flexible member; and 

c) a computer on the vehicle and operatively connected to the vibration sensor and to the sensor which directly 
or indirectly measures the speed, said computer computing the power transmitted by the elongate flexible 
member using the measured vibration and directly or indirectly measured speed of the elongate flexible mem- 
io ber. 

44. The apparatus of claim 43, further including a display for providing a visual display of the computed power. 

45. The apparatus of claim 43, further including signal processing means for eliminating spurious information due to 
is noise sources. 

46. A method for measuring the power transmitted by an elongate flexible member in a system having a first engage- 
ment means separated from a second engagement means, the elongate flexible member coupling the first en- 
gagement means with the second engagement means so that when the first engagement means is rotated by a 

20 power source the elongate flexible member moves in the direction of the rotation of the first engagement means, 

which in turn causes the moving elongate flexible member to rotate the second engagement means; the method 
comprising the steps of: 

measuring directly or indirectly the speed variation of the elongate flexible member; and 
25 receiving the direct or indirect speed variation measurement in a computer, 

generating a signal representing the power transmitted by the elongate flexible member. 

47. The method of claim 46, characterized in that extraneous signals are filtered by analog signal processing. 

30 48. The method of claim 46, characterized in that extraneous signals are filtered by digital signal processing. 

49. The method of claim 46, characterized in that the power output of the elongate flexible member is displayed by a 
digital display of the apparatus. 

35 50. The method of claim 46, characterized in that the signal representing the power output of the elongate flexible 
member is provided by the computer to a external unit coupled to said computer, said external unit displaying the 
power output. 

51. The method of claim 46, characterized in that the system is a bicycle. 

52. The method of claim 46, characterized in that the elongate flexible member is a chain. 

53. The method of claim 46, characterized in that the elongate flexible member is a belt having a plurality of metal 
material equally spaced within and throughout the belt. 

54. The method of claim 46, characterized in that system comprises a rotateable wheel being connected to said the 
second engagement means, said elongate flexible member arranged to rotate said wheel by rotating said second 
engagement means. 

so 55. The method of claim 46, characterized in that the first and second engagement means each comprise various 
sized sprockets so that the elongate flexible member moves from one sprocket to another sprocket on each en- 
gagement means. 

56. The method of claim 46, characterized in that the vibration frequency of the elongate flexible member is measured 
55 to determined the absolute power level at specified time intervals 

57. The method of claim 54, characterized in that the speed variation signal of the elongate flexible member is deter- 
mined by using cadence variation signal or wheel speed variation signal. 
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58. An apparatus for measuring the power transmitted by an elongate flexible member, comprising: 

a system having a first engagement means separated from a second engagement means, the elongate flexible 
member coupling the first engagement means with the second engagement means so that when the first 
s engagement means is rotated by a power source the elongate flexible member moves in the direction of the 

rotation of the first engagement means, which in turn causes the moving elongate flexible member to rotate 
the second engagement means; 

a sensor that directly or indirectly measures the speed variation of the elongate flexible member; and 
a computer that determines the power transmitted by the elongate flexible member by receiving the vibration 
10 measurement and. the speed variation measurement directly or indirectly measured. 

59. The apparatus of claim 58, characterized in that extraneous signals are filtered by analog signal processing. 

60. The apparatus of claim 58, characterized in that extraneous signals are filtered by digital signal processing. 

15 

61. The apparatus of claim 58, characterized in that the apparatus has a digital display for indicating the power output 
of the elongate flexible member. 

62. The apparatus of claim 58, characterized in that the computer generates a signal representing the power output 
20 of the elongate flexible member to an external unit coupled to said computer, said external unit displaying the 

power output. 

63. The apparatus of claim 58, characterized in that the system is a bicycle. 

25 64. The apparatus of claim 58, characterized in that the elongate flexible member is a chain. 

65. The apparatus of claim 58, characterized in that the elongate flexible member is a belt having a plurality of metal 
elements equally spaced within and throughout the belt. 

30 66. The apparatus of claim 58, characterized in that the system comprises a rotateable wheel being connected to said 
the second engagement means, said elongate flexible member arranged to rotate said wheel by rotating said 
second engagement means. 

67. The apparatus of claim 58, characterized in that the first and second engagement means each comprise various 
35 sized sprockets so that the elongate flexible member moves from one sprocket to another sprocket on each en- 
gagement means. 

68. The apparatus of claim 58, characterized in that the vibration frequency of the elongate flexible member is meas- 
ured to determinate the absolute power level at specified time intervals. 



40 



45 



50 



55 



69. The apparatus of claim 58, characterized in that the speed variation signal of the elongate flexible member is 
determined by using cadence variation signal or wheel speed variation signal. 
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variation of the flexible member 



3. Claims: 58-69 

Apparatus for measuring the power transmitted by an elongate 
flexible member comprising a sensor that measures speed 
variation and a vibration measurement 



4. Claims: 34-36 

Method for determining the vibrational frequency of a 
elongate flexible member on a vehicle 
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